Slow convergence tunes onset of strongly discontinuous explosive percolation 
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Contrary to initial beliefs, random graph evolution under an edge competition process with fixed 
choice (an Achlioptas process) seems to lead to a continuous transition in the thermodynamic limit. 
Here we show that a simpler model, which examines a single edge at a time, can lead to a strongly 
discontinuous transition and we derive the underlying mechanism. Starting from a collection of n 
isolated nodes, potential edges chosen uniformly at random from the complete graph are examined 
one at a time while a cap, k, on the maximum allowed component size is enforced. Edges whose 
addition would exceed size k can be simply rejected provided the accepted fraction of edges never 
becomes smaller than a decreasing function, g(k) = 1/2 + (2k)~ 13 . If the rate of decay is sufficiently 
small (/3 < 1), troublesome edges can always be rejected, and the growth in the largest component 
is dominated by an overtaking mechanism leading to a strongly discontinuous transition. If /3 > 1, 
once the largest component reaches size n 1 ^, troublesome edges must often be accepted, leading to 
direct growth dominated by stochastic fluctuations and a "weakly" discontinuous transition. 

PACS numbers: 64.60.ah, 64.60.aq, 89.75.Hc, 02.50.Ey 



Percolation is a theoretical underpinning for analyz- 
ing properties of networks, including epidemic thresholds, 
vulnerability and robustness [TJ-fB] , with large-scale con- 
nectivity typically emerging in a smooth and continuous 
transition. A new class of percolation was conjectured 
to exist two years ago, when it was shown that includ- 
ing a small amount of fixed choice in what to connect up 
leads to an abrupt and seemingly discontinuous "Explo- 
sive Percolation" transition |7J. In just the past months 
it has been established that any random graph evolu- 
tion process with fixed choice in fact leads to a continu- 
ous transition in the thermodynamic limit [81-111]. Here 
we show that a simple stochastic process can lead to a 
discontinuous transition and we analytically derive the 
underlying mechanism. Starting from a collection of iso- 
lated nodes, potential edges chosen uniformly at random 
are examined one at a time while a cap on maximum 
allowed component size is enforced. Undesirable edges 
are simply rejected provided that the accepted fraction 
of edges is never smaller than a decreasing function. We 
show that if this function decreases sufficiently slowly 
then undesirable edges can always be rejected leading to 
growth via overtaking and a strongly discontinuous per- 
colation transition. 

A prototypical process of percolation on a network be- 
gins from a collection of n isolated nodes with edges con- 
necting pairs of nodes sequentially chosen uniformly at 
random and added to the graph [IS]. In a process with 
fixed choice, as studied in [7J, instead, a fixed number 
of randomly selected candidate edges are examined to- 
gether, but only the edge that maximizes or minimizes a 
pre-set criteria is added to the graph. The rule in [7J min- 
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imizes the product of the sizes of the two components to 
be merged. Numerical evidence indicates this rule leads 
to an abrupt jump in the fraction of nodes in the largest 
component at the critical point due to the addition of one 
single edge. Since the publication of [JJ several related 
processes have been shown to result in similar transi- 
tions [T3rH§] . yet, it has also been observed that the size 
of the jump at the critical point decreases with system 
size UHl [SU] (termed "weakly" discontinuous), with 
the most recent results showing that fixed choice per- 
colation transitions on networks are in fact continuous, 
with no jump, in the thermodynamic n — ¥ oo limit [5llll|. 
(Yet, fixed choice processes on lattices may be discontin- 
uous [21].) 

The basic model we analyze here was originally intro- 
duced by Bohman, Frieze and Wormald (BFW) [SSJ , and 
predates the work on explosive percolation. The BFW 
process is initialized with a collection of n isolated nodes 
and a cap on the maximum allowed component size set to 
k = 2. Edges are then sampled one at a time, uniformly 
at random from the complete graph. If an edge would 
lead to formation of a component of size less than or equal 
to k it is accepted. Otherwise the edge is rejected pro- 
vided that the fraction of accepted edges remains greater 
than or equal to a function g(k) = 1/2 + (2k)~ 1 ^ 2 . If the 
accepted fraction would not be sufficiently large, the cap 
is augmented to k + 1 repeatedly until either the edge can 
be accommodated or g(k) decreases sufficiently that the 
edge can be rejected. (For the explicit BFW algorithm 
see the Supplementary information (SI).) 

Here we modify the original function such that g{k) = 
min(l,l/2 + (2k)~@). The parameter /3 controls the rate 
of convergence to the asymptotic limiting value of 1/2, 
as shown in Fig. [lja). Letting d denote the fraction of 
nodes in the ith largest component, we show both analyt- 
ically and via numerical investigation that for f3 < 1 any 
significant growth in C\ is dominated by an overtaking 
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FIG. 1: Convergence rate determines growth mechanism, (a) The function g(k) for different /3. Note, g(k) — 1 means all edges 
must be accepted equivalent to Erdos-Renyi [12] • (b-c) Evolution of Ci and C2 versus edge density, t/n for n = 10 6 . (b) For 
f3 — 1/2, two giant components emerge simultaneously. Inset is the behavior in the critical region showing growth via the 
overtaking process when what was Ci becomes Ci- (c) A typical realization for f3 — 2.0. Inset shows direct growth, that Ci 
and C2 merge together, and what was C3 becomes the new Ci- 



mechanism where smaller components merge together to 
become the new largest component, leading to an evolu- 
tion that converges to a well-defined probability distribu- 
tion and a strongly discontinuous transition. In contrast, 
if > 1 significant direct growth of C\ is allowed and the 
process is dominated by stochastic fluctuations, leading 
to a weakly discontinuous transition that is likely con- 
tinuous as n — > 00. The typical evolution of C\ and C2 
for j3 — 0.5 and (3 = 2.0 are shown in Figs, [ljb) and (c) 
respectively. (The simultaneous emergence of multiple 
stable giant components was recently shown in [23], but 
the underlying mechanism leading to the discontinuous 
transition, our current focus, was not identified.) 

We numerically measure G\ , C% , and C3 throughout 
the evolution process for various (3 € [0.5, 00], for a large 
ensemble of realizations and range of system sizes n. For 
each realization we define the critical point as the single 
edge t c whose addition causes the largest change in the 
value of Ci, with this largest change denoted by AC^ ax . 
As shown in Fig.^a), for f3 — 0.5, AC,j lax is independent 
of system size n and strongly discontinuous. The same 
holds for AC^ lax , the largest jump in Ci- For (3 = 2.0, in 
contrast, the transition is weakly discontinuous meaning 
that AC,^ ax decreases with n. 

As discussed in [5D] , whenever a single edge is added to 
the evolving graph, C\ may increase due to one of three 
mechanisms: (1) Direct growth, when the largest com- 
ponent merges with a smaller one; (2) Doubling, when 
two components both of fractional size G\ merge (this 
is the largest increase possible); (3) Overtaking, when 
two smaller components merge together to become the 
new largest. In [20] it is strictly proven that if direct 
growth is prohibited up to the step when only two com- 
ponents are left in the system, then a strongly discontin- 
uous transition ensues. We next show via analytic argu- 
ments that throughout the subcritical regime, for f3 < 1, 
direct growth only occurs when the largest component 
merges with an essentially isolated node and all signifi- 



cant growth is due to overtaking. In contrast, for j3 > 1 
the evolution is initially the same, but once C\ ~ n 1 ^ , 
large direct growth of C\ dominates. 

Using the notation of [22], let t denote the number of 
accepted edges and u the total number of sampled edges. 
Thus for any k (the maximum allowed component size), 
the fraction of accepted edges t/u > g(k). If t/u is suf- 
ficiently large an edge leading to C\n > k can be simply 
rejected. In contrast, if t/u = g{k) and the next edge 
sampled, denoted e„+i, would lead to G\n > k that edge 
cannot be simply rejected since t/(u + 1) < g{k). One of 
two situations must happen, either: (i) k increases until 
the edge e u+1 is accommodated, or (ii) g{k) decreases suf- 
ficiently that t/(u+ 1) > g(k) and e u +\ is rejected. Thus 
we need to determine the order of the smallest augmen- 
tation of k that makes t/(u + 1) > g(k). For f3 > 0.5 and 
any k the smallest fraction of accepted edges is 

5-«<*>-i + (5)'- (1) 

Rearranging Eq. [T] and differentiating by k yields 

du _ f3 f_ 

dk 20 [1/2 + (l/2fc)' 3 ] 2 k^ 1 

At some point in the subcritical regime we know that 
t ~ 0(n). For /? = 0.5 this has been established rig- 
orously, specifically that by the end of stage k = 25, 
t/n — > 0.841 as n — > 00 [32]. For other /3, we establish 
this via numerical simulation (see SI, Fig S3(a)). With 
this is place we find that once k ~ n 1 (with 7 < 1 for 
the subcritical region) then du/dk ~ n 1_7 ^ _7 . Thus an 
increase in k of order C(n 7/3+7_1 ) is sufficient to ensure 
t/ (u+1) > g(k) and that edge e u+ i can be rejected. But 
there are different behaviors for j3 > 1 and f3 < 1. 

For j3 > 1 there are three regimes. For 7 < l/((3 + 1) 
then 7/3+7— 1 < so the necessary 0(rf / ° +1 ~ 1 ) increase 
in k requires only k k + 1. Then once + 1) <7< 
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FIG. 2: Slow convergence leads to a strongly discontinuous transition and growth by overtaking, (a) For j3 — 0.5, AC^ax and 
A(7„„ are independent of n and both largest components emerge in strongly discontinuous phase transitions. For j3 = 2.0, 



AC' 



showing a weakly discontinuous phase transition, (b) For f3 — 0.5, main plot is T(n) the number of times Ci 



undergoes direct growth versus n, with two regimes separated by k = n 1// ^ +1 '. (Inset) S(n), the average size of component 
that merges with C\ during direct growth, is essentially a constant: S(n) m 1.1 in both regimes, (c) For /3 = 2.0, main plot 
is S(n), showing three regimes (see SI, Fig. SI for /3 — 3). For k < n 1 '^ +1 ' we observe S(n) ~ 1. The intermediate regime 



is noisy. Then once k > n 1 / 8 , random edges must be accepted at times and S(n) 



(Ci merges with other essentially 



macroscopic components). (Inset) T(n) is essentially constant: T(n) « 5. All data points are the average over 30 to 100 
independent realizations (based on system size), with error bars smaller than the symbols unless otherwise indicated. 



1//3, an increase in k of 0(?i 7 ' 3+ ' 7 ~ 1 ) < k ~ C\n is re- 
quired. However, once 7 > 1//3, then C(n 7/3+7_1 ) > 
C\n. The required increase in k is greater than C\n, al- 
lowing C*x to even double in size, and edge e M +i must be 
accepted. So, once in the regime C\n ~ n 1///9 every time 
t/u = g{k) the next edge, e u+ i, must be accepted. In this 
situation, the probability two components are merged be- 
comes, as in Erdos-Renyi [T^], proportional to the prod- 
uct of their sizes and no mechanism remains enforcing 
that the components stay of similar sizes. 

For /3 s [0.5,1) there are only two regimes. Here again 
for 7 < l/(/3 + 1) then k — >• k + 1 allows edge e u+ i to 
be rejected. This regime extends until 7 > l/(/3 + 1), 
when 7/3 + 7— 1 > 0, but now we use the less strict prop- 
erty that /37 + 7 - 1 < /3 and thus n 13 ^^ 1 < < k. 
So k increases in increments of "nP at most, and edge 
e u+ \ is rejected. The slow increase results in multiple 
components of size similar to Ci throughout the evo- 
lution. In particular, for C\ = Sn with 6 <C 1, there 
exist many components of size 0{n). Order them as 
C\n > C2n... > Cin. Assuming " > " strictly holds 
(i.e., choosing non-degenerate components), due to the 
0(nP) increase in k, there will be components such that 
Cj + Cf_i > Cj (or else two smallest components of size 
0(n) would merge very quickly) and there will be a point 
when (Ci + C/_i)n < k < (Cj + C/_2)n, allowing for 
growth by overtaking. We explicitly observe this over- 
taking process for /3 = 0.5 until only three components 
of size 0(n) remain. Then we add edge t c and C2 and 
C3 merge for the final overtaking event as shown inset in 
Fig. [ljb). 

We confirm these predictions via numerical simulations 
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FIG. 3: (a-d) Distribution in values of Ci and C2 at t c ob- 
tained over 100 independent realizations for n = 10 6 . (t c is 
defined as the single edge whose addition causes the biggest 
increase in C\.) (a) Ci for f3 = 0.5. (b) C 2 for j3 = 0.5. (c) 
Ci for /3 = 2.0. (d) C2 for j3 — 2.0. Inset to (c) shows average 
size of Ci and C2 at t c over 100 realizations for different /3. 



using two choices, /3 = 1/2 and f3 = 2.0. Let S(n) denote 
the average size of the component Ci which connects to 
C\ via direct growth for a system of size n, and let T(n) 
denote the number of times direct growth occurs. Fig- 
ure |2fb) is for /3 = 1/2 for the two regimes which are 
separated by k = n 1// ^' 3+1 ^ = n 2 / 3 . As shown in the in- 
set, throughout both regimes S(n) ss 1.1 is essentially a 
positive constant. But T(n) (the main figure) shows a 
distinct regime change. At first T(n) ~ n 019 . Then in 
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FIG. 4: No evident scaling behaviors for /3 = 0.5 whereas quantities for /? = 2.0 exhibit critical scaling, (a) Ci and C2 versus 
t c — t. For ft = 2.0 both Ci, C2 ~ (t c — i) — 3 " 17 ) yet /3 = 0.5 shows no obvious scaling. Inset is the local slope estimate for Ci. (b) 
Distribution of component density n(s) (number of components of size s divided by ri) at different points in the evolution for 
P = 0.5. Inset is n(s) at t/n = 0.96 for various n, showing no finite size effects in the location of the right hump, (c) Evolution 



of n(s) for P = 2.0, with n(s) 



at t c - Inset is n(s) at t/n = 0.93 for various n, again showing no finite size effects. 



the regime starting with k = n 2 / 3, up to and including t c 
we see a much more rapid increase T{n) ~ n . So we 
see direct growth occurring more frequently in the sec- 
ond regime, but the direct growth continues to be due to 
merging with an essentially isolated node. 

Figuregc) is for f3 = 2, with S(n) for the three regimes 
shown in the main plot. Up until k = n}'^ +1 ' = n 1 ^ 
the behavior is the same as for /? = 0.5 as expected since 
Afc = 1 is enough for an edge to be rejected and we see 
S(n) sa 1.1. Then in the regime n 1 / 3 < k < n 1 / 2 , S(n) is 
larger and has large fluctuations. Finally in the regime 
from k = n 1 /* 3 = n 1 / 2 up until edge t c , we see C\ grow in 
large bursts, with S(n) ~ n°", so C\ merges with other 
essentially macroscopic components. (Fig. SI of the SI 
shows equivalent behavior for f3 — 3.0.) As shown inset, 
T(n) is essentially independent of regimes, with T(n) rs 5 
in the first regime and in the third, with T(n) < 1 but 
fluctuating in the second (not shown). 

For p = 2.0, due to the acceptance of random edges 
once k > n 1 / 2 , we observe large fluctuations in the size of 
the giant components at t c . And due to the slow increase 
in k for p — 0.5, the sizes of the components evolve in 
a predictable manner. Figure [3] shows these behaviors, 
with (a) and (b) C\ and Ci observed at t c over 100 in- 
dependent realizations for /3 = 0.5 and (c) and (d) the 
equivalent for ft = 2. Our theoretical analysis holds for 
all P < 1, yet, with numerical analysis we could only con- 
firm these behaviors for /? € [0.5, 0.7], which may be due 
to finite size effects. 

We also carry out a scaling analysis on the general 
BFW model. We first consider the behavior of both C\ 
and Ci near t c , as shown in Fig. Qa). For P = 2 we find 
traditional critical scaling, that C\,Ci ~ (t c — t)~ v with 
7] = 1.17, the same scaling as displayed by the Product 
Rule (PR) edge competition rule studied in [7] [T7J • 
For P = 0.5, Ci, Ci show no obvious scaling behavior. 



More importantly we study the component size den- 
sity n(s) (the number of components of size s divided 
by n). We measure the distribution of n(s) at different 
points in the evolution up to the critical point. Note for 
P = 0.5, t c w 0.976n and for p = 2.0, t c w 0.951n (see 
SI, Fig. S4(a)). For p = 2.0, Fig. gc), the behavior is 
similar to that for PR and other edge competition mod- 
els with fixed choice, where at the critical point there 
is clear scaling behavior, n(s) ~ s~ T with r = 2.1 (the 
same r as for PR). Yet, as shown in Fig. |4]jb), the evo- 
lution for P — 0.5 does not show any scaling. There is 
a pronounced right-hump which forms early in the evo- 
lution, then moves rightward due to overtaking. Inset to 
Fig. gb) and (c) is n(s) at t = 0.96 and t = 0.93 respec- 
tively, for many different n. Note that the peak of the 
right-hump is independent of n. This is quite distinct 
from Erdos-Renyi and PR which show finite size effects; 
for PR the location of the peak moves rightward with 
n, as recently shown in |llj where a finite size scaling 
function for PR is established. 

Also interesting is the rank-size component distribu- 
tion at t c . Ci for p = 0.5 does not show scaling, yet for 
p = 2 we find C t - i~ s with 8 = 0.62 (see SI, Fig.S4(d)). 
(Note the p — 2 model and PR have the same r\ and r 
exponents, but for PR S — 0.90.) Finally we define as 
number of distinct stages visited up to time t c . We find 
that N k - n - 603 for p = 0.5 and N k - n - 292 for p = 2.0 
(see SI, Fig.S4(c)). More stages are visited for p = 0.5 
consistent with our theoretical analysis. 

In summary, we have shown that the convergence rate 
of the generalized BFW model tunes the onset of discon- 
tinuous percolation. For sufficiently slow convergence, 
troublesome edges can always be rejected and significant 
growth in C\ occurs only via overtaking, leading to the 
discontinuous emergence of multiple giant components. 
Algorithms that generate multiple giant components may 
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have a range of applications, such as controlling gel sizes 
during polymerization [25] or creating building blocks for 
modular networks. In addition, the BFW model may be 
simpler to implement in practice than fixed-choice edge 
competition models. A more rigorous understanding of 
the evolution of n(s) for (3 < 1 as n — > oo would be 
useful. Of course real-world networks such as social net- 
works or the Internet, are by definition of finite size, 
so the implications of asymptotic theory may not ap- 
ply. For instance, for the fixed-choice edge competition 
model in [5] it was shown that systems as large as 10 18 
nodes still exhibit a large macroscopic jump of 10% of the 



system size. (PR decays more quickly. Using the scal- 
ing reported in [HI [20l [23] for systems of 10 18 the jump 
would be approximately 3%, for 10 9 it would be 10%.) 
Finally, the Gaussian model 19J which also only controls 
the largest cluster size also appears to be strongly dis- 
continuous |26j . 
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Supplementary Information: 
Slow convergence tunes onset of strongly discontinuous explosive percolation 

Wei Chen, Raissa M. D'Souza 

In the main manuscript we study the BFW model with the generalized acceptance function g(k) = 
min (l, 1/2 + (2fc) _/3 ). We show that if j5 > 1 eventually in the subcritical regime once, k ~ n 1 ^, there are 
times that an undesirable random edge must be accepted. In contrast, if j3 < 1 undesirable edges can always 
be rejected, leading to growth by overtaking and a strongly discontinuous phase transition. Here we present 
additional technical details that are not necessary for the main manuscript but either highlight interesting 
aspects of the BFW model or provide further clarification. 

S.l The BFW algorithm 

Stating the BFW model explicitly, it begins with a collection of n isolated nodes and proceeds in phases 
starting with k = 2. Edges are sampled one at a time, uniformly at random from the complete graph induced 
on the n nodes. Using the notation in [2], let u denote the total number of edges sampled, A the set of 
accepted edges (initially A — 0), and t = \A\ the number of accepted edges. At each step u, the selected 
edge e u is examined via the following algorithm (the BFW algorithm): 

Set I = maximum size component in A U {e u } 
if (/ < k) { 

A <- A U {e„} 

u <- u + 1 } 
else if (t/u<g(k)){ k <- k + 1 } 
else { u «— u + 1 } 

Walking through the loop, if adding e u leads to a component of size I < k, the edge is added and step u 
completes and we move onto step u + 1. Otherwise, check if t/u < g(k). If not, then t/u > g(k) and e u is 
ignored (i.e. rejected) and we move onto step u + 1. If instead it is the case that t/u < g(k) then we augment 
k and walk through from the "if loop from the start again. In other words, while t/u < g(k), k is augmented 
by one repeatedly until either k becomes large enough that edge e u is accepted or g(k) decreases sufficiently 
that edge e u can be rejected at which point step u finally ends. Note g(k) = 1 requires that all edges be 
accepted, equivalent to Erdos-Renyi [1]. Also we use the "min" function, g(k) = min (l, 1/2 + (2fc) - ' 9 ), to 
ensure g(k) < 1 for all choices of /?, since the fraction of nodes accepted cannot exceed unity. 

T. Bohman, A. Frieze, and N. C. Wormald established rigorous results whereby setting p(200) = 1/2, 
all components are no larger than k = 200 nodes (i.e., no giant component exists) when m = 0.96689n 
edges out of 2m sequentially sampled random edges have been added to graph [2]. They further establish 
that a giant component must exist by the time m = c*n out of 2m sampled edges have been added, with 
c* 6 [0.9792,0.9793]. Yet, they did not analyze the nature of the percolation transition. 

In [3] it was shown that the original BFW model leads to the simultaneous emergence of two giant 
components (each with fractional size C > 0.4), and the stability of the multiple giants was analyzed. 
(Essentially, once in the supercritical regime, there are always sufficient edges internal to components sampled 
that whenever an edge connecting two giant components is sampled it can be rejected.) It was also shown 
in [3] that by using the function g(k) = a + (2k)~ 1 / 2 , then the parameter a tunes the number of stable 
giants that emerge simultaneously. 



7 



S.2 Large fluctuations for /3 > 1 

As discussed in the main manuscript, for f3 > 1 once 

k ~ n 1 /' 3 , then there are situations in which an edge must 10 7 

be accepted, regardless of the size of the resulting component. 

(Figure SI shows this for ft = 3.) This leads to large fluctu- 10 

ations in the size of the resulting component, and in the sizes j 

of the components at t c as shown in Fig. 3 (c) and (d) of the 

main manuscript. The typical scenario is as illustrated in Fig. _ 1 ° 4 

1(c) of the main manuscript, that C\ « 0.8 and C*2 emerges in <» . 

a continuous transition starting from Ci = 0. Yet, as seen in 

Fig. 3 (c) and (d), we also observe scenarios where G\ is much i° z 

smaller and Ci is already of macroscopic size at t c . Figure S2(a) , 

shows such a scenario, where C± ~ 0.68 and C2 ~ 0.12 at t c , 

with the inset showing that C^ still emerges via direct growth. 10° , 

10 10 10 10 10 

The situation for p = 1/2, instead is simpler. C\ and Ci have a n 

well defined probability distribution at t c , as shown in the main 

manuscript in Fig. 3 (a) and (b), where C\ = 0.570 ± 0.001 Figure SI: S(n) for (3 = 3, analogous to Fig. 
and C 2 = 0.405 ± 0.001. (Note the error bars were obtained 2(c) in the main, showing d merges with 
in Ref. [3]). Figure 1(b) illustrates this only scenario, with the macroscopic components once k > n 1 ^. 
inset showing that C\ emerges due to overtaking. 

We numerically study the behavior for a few more choices of /3. The inset to Fig. 3 (c) in the main 
manuscript shows the average value and corresponding error bars for C\ and C2 over 100 different network 
realizations at t c for a range of values of [3. Here, Fig. S2(b) shows the corresponding standard deviation in the 
size of these two largest components, denoted as \i > X2 respectively, over a larger ensemble of 800 realizations. 
We also numerically measure the number of times the two largest components ever merge together up to 
and including edge t c for different choices of /3, as show in Fig. S2(c). The interval ft £ (0.7,1) separates 
the region where growth by overtaking dominates from the region where direct growth dominates. Although 
our analysis in the main manuscript holds for f3 < 1, we could only confirm the behaviors numerically for 
P G [0.5, 0.7), which we suspect is due to finite size effect once /3 > 0.7, when occasionally C\ and C2 merge. 



k<n" 4 

n" 4 <k<n 1 ' 3 

k>n 1 ' 3 



[1=3.0 



S(n)~n" 



S{n)~ n 0J 




t/n |i P 



Figure S2: System size here is n = 10 6 . (a) The sizes of C\ and C2 at t c exhibit large fluctuations for /3 = 2.0 
(Fig 3 of main manuscript). Typically C'2 = 0, but there are also significant number of instances where C2 is 
already of macroscopic size at t c , with a typical such realization shown here. Inset shows that, regardless, C\ 
is augmented via direct growth with C2 , and what was the third largest component (which has macroscopic 
size) becomes the new Ci- (b) Standard deviation of C\ (blue line) and C2 (red line) at t c obtained over 800 
independent realizations for each /3, showing a sharp peak at /3 « 1. (c) The total number of times that C\ 
and C2 ever merge together measured over 100 realization, showing that this never happens for (3 < 0.7. 
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Figure S3: (a) Fraction of added edges, t/n, by the end of stage k = 25 versus /3 for different system sizes 
showing a peak at [3 w 1. Values marked n = oo were obtained from linear extrapolation, (b) Fraction 
of added edges, t/n, by the end of stage k = 25 versus different system sizes for ER, PR, and BFW with 
ft = 0.5 and /3 = 2.0. (c) Size of largest component versus different system size at the time when a specified 
fraction of edges have been added for ER, PR, and BFW with /? = 0.5 and j3 = 2.0. 



5.3 Establishing t ~ 0{n) and lack of finite size effects 

To analyze Eqn. 2 in the main manuscript we must establish that t ~ 0(n) in the subcritical regime. 
In [2] they establish this rigorously for /3 = 0.5, specifically that t/n — > 0.841 as n — > oo by the end of 
stage k = 25. Here we establish this fact for a range of f3 and at the same time show the lack of finite size 
effects for the BFW model. Let ^(fc) be the fraction of edges added by the end of stage k. We measure 
-{k = 25) for a series of values of /3 with 0.5 < /3 < 10 for different system sizes, as shown in Fig. S3(a). 
(Each data point is the average over 100 independent realizations.) From these measurements for different 
n we obtain the asymptotic n — > oo value of -(k = 25) by linear extrapolation, showing t/n asymptotically 
converges to positive constant for all ft > 0.5 and the number of added edges t ~ 0(n) at stage k > 25. 
Figure S3(b) shows that = 25) is essentially independent of system size and converges to a positive 
constant for BFW for both f3 = 0.5, 2.0. However, -(fc = 25) decreases to asymptotically for ER and PR, 
with -(fc = 25) ~ n~ T , r = 0.070 and 0.015 respectively. (Note PR is the Product Rule graph evolution 
process studied in [4].) 

On the other hand, rather than measuring t/n for fixed k, we can measure the value of k (which is 
essentially equal to C\) at the time when t/n attains a specified value. Figure S3(c) shows that for BFW 
with /? = 0.5,2.0, k w C\n is essentially a positive constant for t/n = 0.9. Whereas for ER and PR then 
C\n ~ n 6 with 6 = 0.175,0.062 respectively. As also seen in Fig. 4(b) and (c) of the main manuscript, the 
component size density distribution n(s) has a right-side hump for the BFW model (both for f3 = 0.5, 2.0) 
whose peak is independent of system size. In contrast, for PR, the location of the peak (which only exists 
in the subcritical region) increases with system size [5, 6]. 

5.4 Scaling behaviors and bounding t c 

We again focus on BFW model with parameter choices j3 = 0.5 and 2.0. To bound the critical window 
for the transitions we implement the numerical method proposed in [4]. We define to(n) as the maximal 
number of edges for which the largest component has size at most n 7 and t% (n) as the minimal number of 
edges for which the largest component has size larger than An where n is system size, and 7 and A are 
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Figure S4: (a) The critical window as bounded by to/n and t\jn showing convergence in first three significant 
figures to 0.976 and 0.951 for /3 = 0.5 and 2.0 respectively, (b) A/n the fractional number of edges added 
between to and t\ is sublinear in n for both /? = 0.5, 2.0, but decays more rapidly for ft = 2.0 showing a 
tighter scaling window, (c) N^/n, the number of distinct values of k visited up to time t c divided by system 
size n, shows scaling behaviors for both j3 = 0.5 and 2.0, with more distinct stages visited for f3 = 0.5. (d) 



Rank-size component distribution for ft = 0.5, 2.0 at t c . For j3 
with S = 0.62. But no scaling can be discerned for fi = 0.5. 



2.0, fitting for 10 < i < 1000 yields C t 



two parameters. A(7, ^4) = t\(n) — to(n) denotes the number of edges added between these two points. 
We set 7 = 1/2, A = 0.5,0.55 for f3 = 0.5,2.0 respectively. Fig. S4(a) shows that, for the first three 
significant figures, both t /n and t\jn approach the same limiting value limn^oo t c /n = 0.976 for (3 = 0.5 
and limn-^oo t c /n = 0.951 for /3 = 2.0. We also numerically measure that A(l/2, 0.5)/n ~ n -°-' i75 for f} = 0.5 
and A(l/2,0.55)/n ~ n ~ 0A37 for f3 = 2.0, which is shown in Fig. S4(b). The model with /3 = 2.0 results in 
a tighter scaling window because in the subcritical region the direct growth mechanism allows the largest 
component to connect to other large components, thus the largest component can grow very quickly in size. 
The same reason can be applied to explain why this evolution process shows a more delayed onset of giant 
component for /3 = 0.5. Although Figs. S4(a) and (b) provide evidence that the transitions are discontinuous 
for ft = 0.5,2.0, the stronger evidence showing ACj as a function of n in Fig. 2(a) of the main manuscript 
reveals that for /3 = 2.0 the transition is "weakly" discontinuous, yet still shows a significant jump in C\ for 
the values of n attainable in simulation. 

We also measure how the number of stages visited scales with system size n for the BFW model with 
j3 = 0.5 and 2.0 Let Nk be the number of distinct stages visited up to time t c . Fig. S4(c) shows that 
Nk ~ n 603 for p = 0.5 and Nk ~ n 292 for ft = 2.0. This is consistent with our theoretical analysis that 
for ft = 0.5, the stage increases slowly in subcritical region while for /3 = 2.0, direct growth makes the stage 
increases very quickly (see main manuscript). Finally, in the main manuscript, we study n(s) the component 
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size density which starts by emphasizing the abundance of the smallest size components. Also interesting is 
the rank-size component distribution at t c , which instead emphasizes the largest components. As show in 
Fig. S4(d), Ci for j3 = 0.5 does not show scaling, yet for /? = 2 we find Ci ~ i~ s with 5 = 0.62. 

S.5 Stability of multiple giant components 

For any j3 € [0.5, oo] it is always observed that two components of significant size exist at some point 
after the critical point, which remain separate throughout the subsequent evolution. Furthermore we find 
C\ > 0.5 and asymptotically C'2 = l — C\. Similar to the mechanism shown in [3] the stability of simultaneous 
giants is due to the high probability of sampling edges internal to the giant components in the supercritical 
regime. Once C\ > 0.5 (as we determine numerically is the case), the probability of sampling a random 
edge which leads to a component of size larger than stage k is at most P = 2Ci(l — C\) < 1/2. Thus the 
fraction of accepted edges over total sampled edges increases to some positive value larger than 1/2 and k 
stops increasing, meaning to two giant components remain distinct throughout the subsequent evolution. 
For more details see [3]. 

References 

[1] P. Erdos and A. Renyi. Publ. Math. Inst. Hungar. Acad. ScL, 5 (17) (1960). 

[2] T. Bohman, A. Frieze, and N. C. Wormald. Random Structures & Algorithms, 25(4):432-449, (2004). 

[3] W. Chen and R. M. D'Souza. Phys. Rev. Lett. 106, 115701 (2011) 

[4] D. Achlioptas, R. M. D'Souza, and J. Spencer. Science, 323, 1453 (2009). 

[5] R. A. da Costa, S. N. Dorogovtsev, A. V. Goltsev, J. F. F. Mendes. Phys. Rev. Lett. 105, 255701 (2010). 
[6] H. K. Lee, B. J. Kim, H. Park, arXiv: 1103.4439. 



